EPJ manuscript No. 

(will be inserted by the editor) 



m ■ 

§ Self-learning Mutual Selection Model for Weighted Networks 

<n : 

o 



O 

m 

6 
o 

o 

Oh 



> 

o 

(N 

in 
o 

O 
Oh 



X 



Jian-Guo Liu 1 , 3, Yan-Zhong Dang 1 , Wen-Xu Wang 2 , Zhong-Tuo Wang 1 , Tao Zhou 2 , Bing-Hong Wang 2 , Qiang Guo 3 , 
Zhao-Guo Xuan 1 , Shao-Hua Jiang 1 and Ming- Wei Zhao 1 

1 Institute of System Engineering, Dalian University of Technology, Dalian Liaoning, 116023, P R China 

2 Department of Modern Physics, University of Science and Technology of China, Hefei Anhui, 230026, P R China 

3 School of Science, Dalian Nationalities University, Dalian Liaoning, 116600, P R China 

Received: date / Revised version: date 

Abstract. In this paper, we propose a self-learning mutual selection model to characterize weighted evolving 
networks. By introducing the self-learning probability p and the general mutual selection mechanism, 
which is controlled by the parameter m, the model can reproduce scale-free distributions of degree, weight 
and strength, as found in many real systems. The simulation results are consistent with the theoretical 
predictions approximately. Interestingly, we obtain the nontrivial clustering coefficient C and tunable 
degree assortativity r, depending on the parameters m and p. The model can unify the characterization of 
both assortative and disassortative weighted networks. Also, we find that self-learning may contribute to 
the assortative mixing of social networks. 

PACS. 05.65.+b Self-organized systems, 87.23.Ge Dynamics of social systems, 87. 23. Kg Dynamics of 
evolution 



1 Introduction 

In recent years, many empirical findings have triggered 
the devotion of research communities to understand and 
characterize the evolution mechanisms of complex net- 
works including the Internet, the World-Wide Web, the 
scientific collaboration networks and so orL [T| | 2||^||4| | y||p) . 
Many empirical evidences indicate that the networks in 
various fields have some common characteristics. They 
have a small average distance like random graphs, a large 
clustering coefficient and power-law degree distribution ^ 
E], which is called the small- world and scale- free char- 
acteristics. Recent works on the complex networks have 
been driven by the empirical properties of real-world net- 
works and the studies on network dynamics [7|l8|l9|ll()| . 
The first successful attempt to generate networks with 
large clustering coefficient and small average distance is 
that of Watts and Strogatz (WS model) . Another sig- 
nificant model proposed by Barabasi and Albert is called 
scale- free network (BA network) 0. The BA model sug- 
gests that growth and preferential attachment are two 
main self-organization mechanisms of the scale-free net- 
works structure. However, the real systems are far from 
boolean structure. The purely topological characterization 
will miss important attributes often encountered in real 
systems. Most recently, the access to more complete em- 
pirical data allows scientists to study the weight evolution 
of many real systems. This calls for the use of weighted 



network representation. The weighted network is often de- 
noted by a weighted adjacency matrix with element u>ij 
representing the weight on the edge connecting node i 
and j. As a note, this paper will only consider undirected 



graphs where weights are symmetric, i.e. 



The 
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strength Si of node i is usually defined as Si = Y^jerU) Wi i> 
where the sum runs over the neighbor node set r(i). But 
in some cases, the sum can not reflect the node strength 
completely. Take the scientific collaboration networks for 
example, the strength of a scientist include the publica- 
tions collaborated with others and the publications writ- 
ten only by himself or herself. Inspired by this idea, the 
node strength is defined as Si = Y^jerU) w ij + Vii where 
rji is node z's self-attractiveness. As confirmed by the em- 
pirical data, complex networks exhibit power-law distri- 
butions of degree P(k) ~ fc -7 with 2 < 7 < 3 [HUH] and 
weight P(w) ~ w~ e [T2I, as well as strength P(s) ~ s~ Q 
[pi] . The strength usually reveals scale- free property with 
the degree s ~ kP, where f3 > 1 [1^11141115] . Driven by new 
empirical findings, Barrat et al. have presented a simple 
model (BBV for short) to study the dynamical evolution 
of weighted networks JS| ■ But its disassortative property 
can not answer the open question: why social networks are 
different from other disassortative ones? Previous mod- 
els can generate either assortative networks |17U18II19| or 
disassortative ones |12U17lf l8 20 21 . Our work may shed 
some new light to answer the question: is there a generic 
explanation for the difference of assortative and disassor- 
tative networks. 
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Previous network models often adopt the mechanism 
that only newly added nodes could select the pre-existing 
nodes according to the preferential mechanism. However, 
the evolution picture ignores the fact that old nodes will 
choose the young nodes at the same time. Inspired by 
this idea, Wang et al. have presented the mutual selection 
mechanism, which leads to the creation and reinforcement 
of connections [20]. But the model ignored the fact that 
every node would enhance its strength not only by creat- 
ing new links to others, but also could by self-learning. In 
this paper, self-learning means that a node enhances its 
strength only by itself without creating new links to oth- 
ers. Inspired by this idea, we propose a weighted network 
model that considers the topological evolution under the 
general mechanisms of mutual selection and self-learning. 
It can mimic the evolution of many real- world networks. 
Our microscopic mechanisms can well explain the char- 
acteristics of scale-free weighted networks, such as the 
distributions of degree, weight and strength, as well as 
the nonlinear strength-degree correlation, nontrivial clus- 
tering coefficient, assortativity coefficient and hierarchical 
structure that have been empirically observed [TT1I12II221 
23,24,25 . Also, the model appears as a more general one 
that unifies the characterization of both assortative and 
disassortative weighted networks. 



2 Construction of the model 

Our model is defined as following. The model starts from 
Nq = m isolated nodes, each with an initial attractiveness 
So. In this paper, so is set to be 1. 

(i) At each time step, a new node with strength sq and 
degree zero is added in the network; 

(ii) Every node strength of the network would increase by 
1 with the probability p; According to the probability 
(1 — p), each existing node i selects m other existing 
nodes for potential interaction according to the proba- 
bility Equ. 0. Here, the parameter m is the number of 
candidate nodes for creating or strengthening connec- 
tions, p is the probability that a node would enhance 
Vi by 1. 



Hi. 



(i) 



where Si = Yljerli) w v + ^ ^ a P an °^ unlinked nodes 
is mutually selected, then an new connection will be built 
between them. If two connected nodes select each other, 
then their existing connection will be strengthened, i.e., 
their edge weight will be increased by 1. We will see that 
m and p control the evolution of our network. 

The evolution of real-world network can be easily ex- 
plained by our model mechanisms. Take the scientific col- 
laboration networks as an example: the collaboration of 
scientists requires their mutual status and acknowledge- 
ments. A scientists would like to collaborate with others, 
whom have strong scientific potentials and long collabo- 
rating history. On the contrary, he may write paper or 
publications only by himself. When he publishes paper 




Fig. 1. (Color online) Numerical results by choosing p = 0.004. 
The data are averaged over ten independent runs of network 
size N = 7000. (a) Cumulative probability strength distribu- 
tion CP(s) with various m. The results are consistent with a 
power-law distribution CP(s) ~ s a . The inset reports the ob- 
tained values by data fitting (full circles) in comparison with 
the theoretical prediction a — 2 + z/[(l — p) 2 m 2 \ (line), (b) 
Cumulative probability degree distribution CP(k) with vari- 
ous m, which demonstrates that the degree distributions have 
power-law tail, (c) Cumulative probability weight distributions 
with various m, which are consistent with the power-law tail 
P(w) ~ it; 8 , (d) To different m, the average strength Sj of nodes 
with connectivity kt . We observe the nontrivial strength-degree 
correlation s ~ k 13 in the log-log scale. 
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Fig. 2. (Color online) Zipf plot of the degree and node strength 
when p = 0.004. 



as the sole author, his strength also increases, which can 
be reflected by r\. For technological networks with traf- 
fic taking place on them, both the limit of resources and 
the internal demand of traffic increment for keeping the 
normal function of the networks may cause the mutual 
selections. 
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Fig. 3. (Color online) The scale of C(k) and k nn with k for 
various m when p — 0.004. The data are averaged over 10 
independent runs of network size N = 7000. 



Fig. 5. (Color online) The scale of C and r with various m and 
p. The data are averaged over 10 independent runs of network 
size N = 7000. 
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Fig. 4. (Color online) Numerical results by choosing m = 5 
with various p. The data are averaged over ten independent 
runs of network size N = 7000. (a)Cumulative probability 
strength distributions CP (a) with various p. The results are 
consistent with a power-law distribution CP [a) ~ s°. (b) Cu- 
mulative probability degree distributions CP(k) with various 
p, which demonstrate that the degree distributions have power- 
law tail, (c) The clustering coefficient C(k) depending on con- 
nectivity k for various p. (d) Average connectivity k nn of the 
nearest neighbors of a node depending on its connectivity k for 
different p. 

3 characteristics of the model 

Considering the rule that my is updated only if node i and 
j select each other, and using the continuous approxima- 
tion, then the time evolution of weight can be computed 
analytically as 



dwij 



(l—p)msj (l—p)msi 



(1 -p) 2 m 2 (sj5j) 



(2) 



Hence, the strength Si(t) is updated by this rate 

(1 — p) 2 m 2 Si 



dsj 
dt 

Notice that 



dt 



1 y ds t 

^\ l dt = 
dt ., 



"(l-p) 2 m 2 s 



P- 



Thus, Equ. © can be expressed by 
(1 — p) 2 m 2 Si 



dsj 
dt 



(1 — p) 2 m 2 t + pt 2 



(3) 



dt. 



(4) 



When p ~ 0(N 1 ), the solution of Equ. can be ob- 
tained approximately as follows 



where 



A = 



Si(t) ~t A , 
(1 - p) 2 m 2 



(5) 



(1 — p) 2 m 2 + z ' 

and z = pN is a constant. Then, we can get that the 
strength distribution obeys the power-law P(s) ~ s~ a 
jlfij with exponent 



a = 1 + i = 2 
A 



(1 — p) 2 w? 



(6) 



One can also obtain the evolution behavior of the weight 
distribution P(w) ~ w e [28], where 



= 2 



2z 



(1 — p) 2 m 2 — z ' 

and the degree distribution P(fc) ~ fc~ 7 , where 

z 



7 



(1 — p) 2 m 2 



(7) 



(8) 



By choosing different values of p and to, we perform nu- 
merical simulations of networks which is consistent with 
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the theoretical predictions. Fig. l(a)-(d) present the prob- 
ability distributions of strength, degree and weight, as well 
as the strength-degree correlation, fixed p = 0.004 and 
tuned by to. Fig. 1(a) gives the node strength distribution 
P(s) ~ s a , which is in good agreement with the theoret- 
ical expression Equ. JBJ. Fig. 1(b) gives the node degree 
distribution P(k) ~ fc -7 . Fig. 1(c) reports the probability 
weight distribution, which also shows the power-law be- 
havior P(w) ~ w e . Fig. 1(d) reports the average strength 
of nodes with degree fej, which displays a nontrivial power- 
law behavior s ~ k@ as confirmed by empirical measure- 
ments ^21- Fig. 2(a)-(b) show the Zipf plot of the simu- 
lation results by fixing a moderate value p = 0.004 and 
varying m. Fig. 2(a) confirms with the math collabora- 
tion network and the Zipf plot of Fig. 1(a) in Rcf. |24| . 
Fig. 3(a)-(b) give the clustering coefficient C(k) depend- 
ing on connectivity k and the average connectivity k nn of 
the nearest neighbors of a node for various to. From the 
numerical results, we can obtain the conclusion that the 
larger the probability p, the larger the effect of exponen- 
tial correction at the head. However, the power-law tail 
which again recovers the theoretical exponent expressions 
can still be observed. Fig. 4 gives the numerical results for 
various p when to = 5. 

Depending on the parameters p and to, the unweighted 
clustering coefficient C, which describes the statistic den- 
sity of connected triples, and degree assortativity r [23 
|2?| are demonstrated in Fig. 5. The assortative coefficient 
r can be calculated from 

where ji , fe« are the degrees of the vertices at the ends of 
the zth edge, with i = 1, 2, • • • , M. From Fig. 5(a), we can 
find that C for fixed m increases with p slightly, and C 
for fixed p monotonously increases with to. The cluster- 
ing coefficient of our model is tunable in a broad range 
by varying both m and p, which makes it more power- 
ful in modelling real-world networks. As presented in Fig. 
5(b), degree assortativity r for fixed p decreased with m, 
unlike the clustering case. While r for fixed m increases 
with p slightly. The model can generates disassortative 
networks for small to and large p, which can best mimic 
technological networks. At large p and small m, assorta- 
tive networks emerge and can be used to mimic social net- 
works, such as the scientific collaboration networks. In the 
model, enhancing the probability p can be considered as 
the probability that a node would like to study by itself to 
enhance its attractiveness or prestige. In the competitive 
social networks, all nodes face many competitors. In order 
to subsist or gain honorableness, they must enhance their 
attractiveness or ability by studying themselves or collab- 
orating with others. This explains the origin of assortative 
mixing in our model and may shed light on the open ques- 
tion: why social networks are different from other disas- 
sortative ones? For example, in the scientific collaboration 
networks, the attractiveness of a scientist could not be 
represented simple by the publications collaborated with 
others. Indeed, there are many other important qualities 



that will contribute to the attractiveness of a scientist, 
for instance, the publications written by himself, etc. Per- 
haps the different self-learning probability contributes to 
human beings fundamental differences. On the other hand, 
to indicates the interaction frequency among the network 
internal components. If m increases, the hubs would link 
more and more "young" sites. Thus, the reason why the 
disassortativity of the model is more sensitive to to lies in 
that collaboration is more important than self-learning in 
the technological networks. In addition, the components of 
technological networks are usually physical devices, which 
can not study by itself. Combining these two parameters 
together, two competitive ingredients, which may be re- 
sponsible for the mixing difference in real complex net- 
works, are integrated in our model. 

4 Conclusion and Discussion 

In summary, integrating the mutual selection mechanism 
between nodes and the self-learning mechanism, our net- 
work model provides a wide variety of scale-free behav- 
iors, tunable clustering coefficient, and nontrivial degree- 
degree and strength-degree correlations, just depending 
on the probability of self-learning p and the parameter m 
which governs the total weight growth. All the statistic 
properties of our model are found to be supported by var- 
ious empirical data. Interestingly and specially, studying 
the degree-dependent average clustering coefficient C{k) 
and the degree-dependent average nearestneighbors de- 
gree k nn {k) also provides us with a better description of 
the hierarchies and organizational architecture of weighted 
networks. Our model may be beneficial for future under- 
standing or characterizing real networks. 

Due to the apparent simplicity of our model and the 
variety of tunable results, we believe our present model, for 
all practical purposes, might be applied in future weighted 
network research. 

This work has been supported by the Chinese Natural 
Science Foundation of China under Grant Nos. 70431001, 
70271046 and 70471033. 
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